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ABSTRACT 



This work investigates some aspects of the transport of 
low-energy e+ (i?k 1 MeV) in the interstellar medium 
(ISM). We consider resonance interactions with magneto- 
hydrodynamic (MHD) waves above the resonance thresh- 
old. Below the threshold, collisions take over and deflect 
e+ in their motion parallel to magnetic-field lines. Using 
Monte-Carlo simulations, we model the propagation and 
energy losses of c+ in the different phases of the ISM 
until they annihilate. We suggest that e+ produced in 
the disk by an old population of stars (SNIa, LMXB), 
with initial kinetic energies below 1 MeV, and propagat- 
ing in the spiral magnetic field of the disk, can probably 
not penetrate the Galactic bulge. 

Key words: Positrons; Diffusion process; Interstellar 
medium. 



1. INTRODUCTION 



The map of the 511 keV e~e"'" annihilation emission, 
measured by the spectrometer SPI on board ESA's INTE- 
GRAL observatory, indicates that c+ annihilate mainly 
in the Galactic bulge, with a bulge-to-disk flux ratio 
B/D ~ 2±l iU. Under the assumption that e+ an- 
nihilate near their sources, this ratio provides severe con- 
straints on the e+ sources. The measured value of B/D 
does not correspond to any known distribution of astro- 
physical objects. The annihilation emission from the disk 
may be due to Al and ''^Ti ejected by massive stars, but 
the origin of e+ in the bulge is still a mystery. Jean et al. 
lITill suggested that e+ produced in the bulge with ener- 
gies lower than a few MeV do not escape and, therefore, 
annihilate in the bulge. Prantzos [16] argued that a bal- 
listic transport of 1 MeV e+ from the disk into the bulge 
via a regular dipolar magnetic field could explain a high 
B/D ratio. 

The transport of c+ is still poorly understood, but, as all 
charged particles, e+ are sensitive to the Galactic electro- 
magnetic field, the flows of interstellar matter and the in- 
teractions with gas particles. In this paper, we report on 



preliminary results of transport in the Galaxy. Sec- 
tion |2] describes the effects of interactions with Alfven 
waves. Section[3]presents the model of coUisional diffu- 
sion. The implications of our model and the other possi- 
ble processes of diffusion are discussed in Section|4l 



2. INTERACTIONS WITH ALFVEN WAVES 



In this section, we restrict the present analysis to the inter- 
actions between c+ and Alfven waves, usually assumed 
to be the main transport agent in the ISM. The effects of 
e"*" interactions with magnetosonic waves are postponed 
to a future work. 

Regardless of the origin of Alfven waves, the strong 
damping of Alfven waves above the proton cyclotron fre- 
quency implies an energy threshold (^^ql) to the reso- 
nance condition. In the case of e+, this condition can be 
satisfied only if Hi" 



7/3 > 12.85 X 10'^ - 

\/'^cm~^ 



(1) 



where 7 is the e+ Lorentz factor, (3 the g+ velocity di- 
vided by the speed of light, B the magnetic-field strength 
and n the total particle density of the considered phase. 
When e+ are in resonance with Alfven waves, we use 
quasilinear theory to derive the diffusion coefficient, 
which enables us to estimate the distance travelled by e+ 
(see Equation 7 in lITill '). 

Note, however, that when the amplitudes of magnetic 
fluctuations become comparable to the strength of the 
mean magnetic field, quasilinear theory gives a crude de- 
scription of e+ transport. This particular aspect will be 
addressed in a future paper. 



3. INTERACTIONS WITH GAS PARTICLES 



Below the quasilinear energy threshold, the resonance 
condition is not satisfied. Therefore, we assume that 
e"*" propagate in a collisional regime. We calculate the 



distance travelled by c+ using Monte-Carlo simulations 
which take into account the physics of c+ interactions 
with gas particles and their associated cross sections 
lEHlHlllH^l- The present analysis is restricted to a neu- 
tral hydrogen gas (n = nni) and a fully ionized hydrogen 
gas (n = n^- + nnii)- The effects of the ionization frac- 
tion and metallicity will be presented elsewhere. Since 
the n-dependence of the quantities of interest is known 
analytically, we may arbitrarily adopt n — 1 cm""^. 

The medium is assumed to be pervaded by a homoge- 
neous magnetic field B^, directed along the z axis and 
having a strength of 5 (in agreement with the typi- 
cal value measured in the Galactic disk). Each e+ travels 
along magnetic-field lines and describes a helical orbit of 
radius equal to its Larmor radius. When a c"*" interacts 
with a gas particle, its trajectory is deflected and the e+ 
loses a fraction of its energy or it annihilates. Each sim- 
ulation is performed using a total number of c+ ranging 
between 5 x 10"^ and 2 x 10"*. The initial pitch angle is 
chosen randomly according to an isotropic distribution in 
the half-space z > 0. In this regime, the distance trav- 
elled along magnetic-field lines by c+ with a given initial 
energy does not depend on the magnetic-field strength. In 
the next sections, the results of our simulations are pre- 
sented for e"*" with an initial kinetic energy of 3 keV. 



3.1. Transport in a neutral medium 

In a neutral medium, c"*" lose energy through ionization 
and excitation of atoms and through elastic scattering 
with atoms. They annihilate with bound either di- 
rectly or after forming a positronium in flight by charge 
exchange. To estimate the distance travelled by e+ in a 
neutral medium, we distinguish three periods during the 
e+ lifetime: 



1. The slowing-down period. 

When the kinetic energy of e"*" lies between the 
quasilinear threshold and the charge exchange 
threshold (6.8 eV in H), e"*" lose energy mainly 
through ionization and/or excitation, and ~ 98 % 
of them annihilate in flight, in quite good agree- 
ment with til]. In Figure [U we show the final dis- 
tribution of e''" parallel to magnetic-field lines. This 
distribution is obtained when the e+ energy reaches 
the charge exchange threshold or when c+ annihi- 
late. The distribution provides an estimate of the 
average distance travelled by e+ along magnetic- 
field lines. It can be fitted by a Gaussian, using 
minimization, with an average distance (z) ~ 
(0.886 ± 0.005)/ncin-3 pc and a spatial dispersion 
(z2)i/2 ^ (0.707±0.003)/ncm-3pc, obtained after 
an average slowing-down time {t) ~ %&/n^^-^-3 yr. 

The distance travelled by c+ during this period does 
not depend on the temperature of the gas, since their 
kinematic energy is always larger than the thermal 
energy of the gas (< 1 eV). 
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Figure 1. Spatial distribution ofe^ parallel to magnetic- 
field lines at the final stage of the slowing-down period. 
The initial e+ kinetic energy is 3 keV and the total gas 
density is n ~ 1 cm^^. 

Note that, in a neutral medium, Alfven waves could 
be strongly damped by ambipolar diffusion. In this 
case, the slowing-down period of collisional diffu- 
sion would start at higher energies, maybe even at 
the injection energy of e+. Starting collisional diffu- 
sion with an initial c+ kinetic energy of 1 McV and 
using the theory of inelastic scattering [6], we obtain 
(z) ~ 4.0/ncin-3 pc, (z^)^/^ ~ 7.4/ncin-3 pc and 
an averaged slowing-down time ~ 22Q / yr. 

2. The thermalization period. 

e"*" below the energy threshold for charge exchange 
lose energy elastically until their kinetic energy 
drops to that of gas particles. We obtain the spa- 
tial distribution parallel to magnetic-field lines that 
is shown in Figure |2] The thermalization time is 
~ lO^/ricni-a yr- Using minimization, we can 
fit the curve by a Gaussian with an average distance 
(z) ~ (4 ± 6) X 10~^/ncni-3 pc and a spatial dis- 
persion (z2)i/2 ^ (0.19±0.01)/n„„-3 pc. During 
this period, c+ annihilation is possible, but negligi- 
ble (< 1%). The effect of the gas temperature is also 
found to be negligible. Between 8000 K and 10 K, 

1/2 

the difference in (z^) is less than 8 %. 

3. The thermalized period. 

During this period, the kinetic energy of e+ is com- 
parable to the thermal energy of gas particles and e"*" 
scatter elastically with atoms. On average, there is 
neither gain nor loss in energy and e+ diffuse over a 
distance A = \/2Dt along magnetic-field lines un- 
til they annihilate directly with bound (t is the 
c+ lifetime |7]). Given in Table[T]are the values of 
the diffusion coefficient D parallel to magnetic-field 
lines that we computed for the different ISM phases. 
It has to be noted that both D and r depend on the 
gas temperature. 

In agreement with the equation of motion of e+ in a mag- 
netic field, the magnetic-field strength influences only the 
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Figure 2. Spatial distribution ofe^ parallel to magnetic- 
field lines at the final stage of the thermalization pe- 
riod. The initial kinetic energy is equal to the charge- 
exchange threshold and the total gas density is n — 
1 cni~^. 
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Table 1. Estimated values of the parallel diffusion co- 
efficient in the thermalized period in the different ISM 
phases. 



perpendicular transport. In all three periods listed above, 
distances traveled perpendicular to magnetic-field lines 
are found to be negligible (d_L < lO^^^/yicm-a pc). 



3.2. Transport in an ionized medium 

In an ionized medium, e+ lose energy through Coulomb 
scattering with ions and free e^. They annihilate with 
free either directly or after forming a positronium by 
radiative combination ItII . Here, we model Coulomb in- 
teractions according to ill ill and 12] • In contrast to propa- 
gation in a neutral medium, propagation here depends on 
the gas temperature. To estimate e+ transport in an ion- 
ized medium, we divide the e+ lifetime into two periods: 




Figure 3. Spatial distribution ofe^ parallel to magnetic- 
field lines at the end of the slowing-down period, in the 
warm ionized ( upper blue line) and the hot ionized (lower 
red line) phases. The initial e+ kinetic energy is 3 keV 
and the total gas density is n — 1 cm~'^ 
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Table 2. Average distance (z) travelled by during the 
slowing-down period and spatial dispersion of the dis- 
tribution the warm ionized and hot ionized 
phases, n is the total gas density. 



The thermalized period. 

During this period, the kinetic energy of e+ is com- 
parable to the thermal energy of gas particles and 
e"*" interact with via Coulomb scattering. On av- 
erage, there is neither gain nor loss in energy and c+ 
diffuse over a distance A = \/2Dt along magnetic- 
field lines until they annihilate with free c^. Here, 
too, the values of the diffusion coefficient D parallel 
to magnetic-field lines are given in Table [T] 



As in the case of a neutral medium, diffusion perpen- 
dicular to magnetic-field lines is found to be negligible 

(di < 2 X 10-iVn„„-3 pc). 



4. DISCUSSION AND CONCLUSION 



1 . The slowing-down or thermalization period. 

During slowing-down, e+ lose energy through 
Coulomb scattering, and a negligible fraction (< 
1%) of e+ annihilate. In Figure[3] we show the spa- 
tial distribution of e+ parallel to magnetic-field lines 
at the end of the slowing-down period, i.e at a time 
^ 25/7icm-3 yi'- The average distance (z) travelled 
by e"*" and the spatial dispersion (z^)^/^ of the dis- 
tribution are indicated in Table |2] for both the warm 
ionized and the hot ionized phases. 



In this section, we discuss the effects of the transport of 
e"*" produced around 1 McV by an old stellar population 
(type la supernovae, low-mass X-ray binaries, ... ) in 
the Galactic disk. Table [3] gives the average distances 
travelled by e+ at different locations in the Galactic disk 
(r > 3 kpc), as computed with the method presented in 
the previous sections. Because e+ follow magnetic-field 
lines, we start by describing the magnetic-field properties 
in the disk. 

The current observational status, based on studies of the 
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Table 3. Average distances travelled by e+ with an initial kinetic energy ofl MeV in the quasi-linear {d^^ and collisional 
{d^^^ ) regimes, at different locations in the Galactic disk. The total particle densities are derived from 0/ 



Galactic synchrotron emission 1 1 ] and of Faraday rotation 
measures |9, 10, 12, 18, 19], is basically the following: 
The interstellar magnetic field has uniform (large-scale) 
and turbulent (small-scale) components of comparable 
strengths (~ 1.5 fiG and ^ 5 fiG, respectively in the 
Solar vicinity). The direction of the uniform magnetic 
field is nearly horizontal in most of the Galactic disk, 
and on the whole, the azimuthal component dominates. 
The uniform magnetic field probably has a spiral pat- 
tern, but it is not known whether this spiral is axisym- 
metric, bisymmetic, or a mixture of different azimuthal 
modes. The strength of the uniform magnetic field in- 
creases smoothly toward the Galactic center, reaching at 
least 4 pG at r ~ 4 kpc. Along the vertical, the uniform 
field strength decreases away from the midplane, proba- 
bly following a two-layer structure with respective scale 
heights ^ 200 pc and ~ 2 - 4 kpc. 

The old stellar population is more spread out along the 
vertical than neutral gas. Here, we assume for simplic- 
ity that at 2 ^ 0.3 kpc the ISM is fully ionized. In 
accordance with the ISM distribution 0, Ht], we adopt 
n ~ 0.2 cm~'^ for the warm ionized phase and n ^ 
0.002 cm~^ for the hot ionized phase. With an initial 
kinetic energy ^ 1 MeV, e+ start propagating by diffu- 
sion on Alfven waves until their kinetic energy falls be- 
low ^5.2 keV in the warm ionized phase and ^ 380 keV 
in the hot ionized phase. While e+ are transported by 
Alfven waves, they travel an average distance ~ 40 pc 
in the warm ionized phase and ^ 250 pc in the hot ion- 
ized phase. Below the resonance threshold, e+ propagate 
by collisions until they annihilate. The average distance 
travelled in the collisional diffusion regime is < 1 pc in 
the warm ionized phase and ^ 600 pc in the hot ion- 
ized phase (see Table O. Consequently, the maximum 
average distance travelled by e+ before annihilation is 
^ 850 pc. We also estimate that the maximum possi- 
ble distance travelled by a typical e"*" is ~ 1 kpc. As- 
suming a spiral magnetic field in the Galactic disk, we 
then conclude that c+ produced at Galactocentric radius 
r > 3 kpc can not reach the Galactic bulge. 

Note that in the current model, e+ follow the average 
magnetic-field lines. However, in the Galactic disk, su- 
perbubbles and supernova remnants generate turbulence 
that can drive magnetic-field lines at high altitude (> 
1 kpc). There, c+ are in a low-density medium where 
complex magnetic effects could significantly change the 
transport of e+. Galactic winds could also transport 6+ 
advectively. In the Galactic center, a number of young 
massive stellar clusters blowing powerful winds have 



been observed ITsl [TtIi . Such points will be addressed 
in a future work. 
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